


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


2005-06 


Detection of frequency-hopped signals 
embedded in interference waveforms 


Brown, Christopher K. 


Monterey California. Naval Postgraduate School 


http://ndl.handle.net/10945/2180 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


/ (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

; | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





NAVAL 
POSTGRADUATE 
SCHOOL 


MONTEREY, CALIFORNIA 


THESIS 


DETECTION OF FREQUENCY-HOPPED SIGNALS 
EMBEDDED IN INTERFERENCE WAVEFORMS 


by 
Christopher Brown 


June 2005 


Thesis Advisor: R. Clark Robertson 
Second Reader: Frank Kragh 





Approved for public release; distribution is unlimited 


THIS PAGE INTENTIONALLY LEFT BLANK 


| REPORT DOCUMENTATION PAGE Form Approved OMB No. 0704-0188 


Public reporting burden for this collection of information is estimated to average | hour per response, including 
the time for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and 
completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington 


headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 


1. AGENCY USE ONLY (Leave blank) | 2. REPORT DATE 3. REPORT TYPE AND DATES COVERED 
June 2005 Master’s Thesis 


4. TITLE AND SUBTITLE: Title (Mix case letters) 5. FUNDING NUMBERS 
Detection Of Frequency-Hopped Signals Embedded In Interference Waveforms 


6. AUTHOR(S) 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING 
Naval Postgraduate School ORGANIZATION REPORT 
Monterey, CA 93943-5000 NUMBER 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING/MONITORING 
N/A AGENCY REPORT NUMBER 





11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 
policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 
Approved for public release; distribution is unlimited. 
13. ABSTRACT (maximum 200 words) 

Many military communications systems utilize frequency-hopped spread spectrum waveforms to protect 
against jamming and enemy detection. These waveforms may be subjected not only to intentional jamming but 
may also be unintentionally jammed by other communications signals. While some systems can overcome in- 
band interference with more signal power, covert systems may be limited to small amounts of transmitted power. 
The objective of this thesis was to investigate a method for resolving a frequency-hopped signal embedded in in- 
terference waveforms. 

With exponential averaging in the frequency domain, the spectra of the interfering signals can be estimated as 
long as they are present over a period longer than that of the frequency-hopped signal. Certain FFT sizes and 


weights are more beneficial to achieving this estimate than others. The interference estimate can be used to extract 
the desired frequency-hopped signal through spectral division of the received signal with the estimate. This tech- 
nique is designated as noise-normalization. Simulations in MATLAB demonstrate the use of the technique and 
show how the desired signal can be resolved. 


14. SUBJECT TERMS 15. NUMBER OF 
Noise-normalization, frequency-hopped, interference estimate, exponential average PAGES 
a7 


16. PRICE CODE 


17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION 
CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF OF ABSTRACT 
REPORT PAGE ABSTRACT 

Unclassified Unclassified Unclassified UL 


NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 





THIS PAGE INTENTIONALLY LEFT BLANK 


il 


Approved for public release; distribution is unlimited 


DETECTION OF FREQUENCY-HOPPED SIGNALS EMBEDDED IN 
INTERFERENCE WAVEFORMS 


Christopher K. Brown 
Ensign, United States Naval Reserve 
B.S., United States Naval Academy, 2004 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN ELECTRICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 


June 2005 
Author: Christopher Brown 
Approved by: R. Clark Robertson 
Thesis Advisor 
Frank Kragh 


Second Reader 


John P. Powers 
Chairman, Department of Electrical and Computer Engineering 


il 


THIS PAGE INTENTIONALLY LEFT BLANK 


1V 


ABSTRACT 


Many military communications systems utilize frequency-hopped spread spec- 
trum waveforms to protect against jamming and enemy detection. These waveforms may 
be subjected not only to intentional jamming but may also be unintentionally jammed by 
other communications signals. While some systems can overcome in-band interference 
with more signal power, covert systems may be limited to small amounts of transmitted 
power. The objective of this thesis was to investigate a method for resolving a fre- 


quency-hopped signal embedded in interference waveforms. 


With exponential averaging in the frequency domain, the spectra of the interfering 
signals can be estimated as long as they are present over a period longer than that of the 
frequency-hopped signal. Certain FFT sizes and weights are more beneficial to achieving 
this estimate than others. The interference estimate can be used to extract the desired fre- 
quency-hopped signal through spectral division of the received signal with the estimate. 
This technique is designated as noise-normalization. Simulations in MATLAB demon- 


strate the use of the technique and show how the desired signal can be resolved. 
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EXECUTIVE SUMMARY 


This thesis investigated recovery techniques for frequency-hopped spread spec- 
trum signals embedded in interference waveforms. Many military communications sys- 
tems utilize frequency-hopped spread spectrum techniques to protect against jamming 
and enemy detection. These systems operate in environments that may contain other 
communications signals and, as a result, may be subjected to interfering waveforms. 
While some systems can overcome in-band interference with more signal power, covert 
systems may be limited to small amounts of transmitted power. The objective of this the- 
sis was to investigate a method for resolving a frequency-hopped signal amongst interfer- 


ence waveforms that might be operating in the same band. 


The spectra of the interfering signals can be estimated with an exponential averag- 
ing algorithm in the frequency domain. Several N-sample windows of the received signal 
are transformed to the frequency domain by taking the fast Fourier transform (FFT). 
These windows of the received signal are averaged to create an N-point frequency do- 
main representation of the interference waveforms. The estimate 1s formed under the 
assumption that the interference waveforms are constant over multiple hops. As a result, 
an individual frequency hop present in an FFT window will not be present in the final 


average of all the FFTs, leaving just the interference waveforms. 


Certain FFT window sizes and weight factors are more beneficial to achieving 
this estimate than others. As the FFT window size increases, the number of FFTs used to 
create the estimate decreases. Too large of an FFT window causes the estimate to be 
formed from a small number of windows, and thus the frequency-hopped signal is not 
averaged out. Small FFT windows cause resolution problems though, and do not accu- 
rately represent the overall shape of frequency spikes and spectral lobes. A balance must 
be found in choosing the FFT window size. Moreover, the weight factor affects how 
much influence an individual window has on the estimate. A large weight value places 
more emphasis on the most recent frequency hop, while a smaller weight accentuates the 
first frequency hop in the sequence. The goal 1s to find a weight that allows each fre- 


quency hop to be reduced to an acceptable level. A good weight factor forms an estimate 


XV 


that 1s an accurate representation of the interference waveforms, with few frequency hops 


remaining in the spectrum. 


The estimate can be used to extract the desired frequency hops through spectral 
division of the received signal with the estimate. This technique is known as noise- 
normalization. Simulations in MATLAB demonstrate the use of the technique and show 
how the frequency-hopped spectra and desired signal can be resolved. A frame-by-frame 
division of the received signal by the estimate reveals frequency hops that are eclipsed by 
the spectral lobe. Given enough frames to create the estimate, a receiver 1s able to re- 
cover all but the last two hops that occur over a finite period. This technique also proves 
to be superior when compared with spectral subtraction. Specific FFT window sizes and 
weight factors lead to clearer normalized spectra and recover more frequency hops than 
other combinations of FFT size and weight factor. The weight factor that provides the 
best estimate is not necessarily the best weight for resolving the frequency-hopped signal. 
Bit error rate simulations of the normalization technique reveal that the noise- 
normalization technique achieves the same bit error rate as a 50-dB reduction in noise 


power. 


The technique described in this thesis, though implemented only in simulation, 
forms an accurate estimate of interference waveforms and successfully detects frequency 


hops eclipsed by other signals in the band. 


XV1 


I. INTRODUCTION 


A. BACKGROUND 

Within military communications, frequency-hopped systems play a major role. 
These systems use changing carrier frequencies so that the signal 1s difficult to intercept. 
Frequency-hopped signals operate within a band, and each adjacent carrier is separated 
by some bandwidth, Af. Given that in some areas the electromagnetic spectrum is not 
strictly regulated, it 1s possible for frequency-hopped radio links to operate amidst other 
signals. In the case of covert transmission, simply increasing transmission power to 
overcome interference that happens to be jamming the receiver 1s not always an option. 
A receiver that has a sense of the environment it is operating in should be able to adjust 


for the interference and recover the desired signal. 


Previous work in [1]-[4] has yielded several different methods for frequency- 
hopped signal detection. These methods all assume that the frequency-hopped signal 1s 
being received in additive white Gaussian noise and thus has a standard mathematical 
estimate. From [5], it is known that the power spectral density of such noise is constant 
across all frequencies. The effect of real communications signals such as phase-shift 
keyed (PSK), minimum-shift keyed (MSK), and quadrature amplitude modulated (QAM) 
signals that might act as interference are more difficult to model. These signals’ spectra 
contain spectral spikes, lobes, and nulls that cannot be approximated as additive white 
Gaussian noise (AWGN). This thesis investigates the detection of a frequency-hopped 
waveform embedded in an interference waveform consisting of a combination of a PSK 
signal and a number of continuous wave (CW) signals. 

B. OBJECTIVES 

The focus of this thesis was two-fold. The first objective was to develop a 
method to estimate the interference environment. The receiver must be able to use this 
estimate for further signal processing. With an N-point window, many small fast Fourier 
transforms (FFT) of the received signal can be averaged to form an estimate of the inter- 
ference power spectral density. The parameters surrounding this method, including the 
size of the FFT window and the configuration of the averaging algorithm, are the key 


components in deriving the interference estimate. 
l 


The second goal was to explore methods to resolve the frequency-hopped wave- 
form using the interference estimate. The goal was not to amplify the desired signal but 
to reduce the noise power level, thus increasing the signal-to-noise ratio. Detection of the 
frequency-hopped waveform was the end goal but was done while keeping in mind that 
the signal contains modulated data that eventually must be recovered. Spectral subtrac- 
tion and noise-normalization were studied. Although noise-normalization has been ana- 
lyzed from a theoretical standpoint, to the best of the author’s knowledge, this 1s the first 
time noise-normalization has been implemented. 

C. ORGANIZATION OF THESIS 

This thesis consists of four parts. Background information on the signals used 
and the MATLAB program utilized for simulating the signal processing are described in 
Chapter II. Chapter [III presents a description of the investigation into developing a tech- 
nique for estimating the interference in the channel. Chapter IV contains a description of 
frequency hop recovery techniques. Chapter V contains a bit error rate simulation for a 
transmitted signal with the recovery technique performed on the channel. The final chap- 


ter, Chapter VI, contains conclusions and recommendations for further research. 


Hl. SIMULATION SIGNALS 


The signals used in detection simulations are described in this chapter. In order to 
test the various aspects of forming an interference estimate and recovering a frequency- 
hopped signal, discrete signal files were used within MATLAB simulations. The 
MATLAB program can import various types of discrete signal files into an array, includ- 
ing files denoted as signals with a .sig extention, which were used for this thesis. 
Mathematical computations, loops, if/else statements, fast Fourier transforms, and bit er- 
ror rate calculations are some of the unique capabilities that the MATLAB program pro- 
vided for creating various simulations. All plots shown in this paper were created using 


arrays in MATLAB. 


The signals used for the stmulation were each sampled at 50 MHz and contain 
1,032,192 points of data. Three files were imported into MATLAB but were manipulated 
in the program to create several different signals for the simulation. Two signals were 
used to create an interference signal that operates in the same band as a transmitted sig- 
nal. The terms “transmitted signal” and “desired signal” refer to the same frequency- 
hopped signal and are used interchangeably throughout this thesis. 

A. DESIRED FREQUENCY-HOPPED SIGNAL 

The desired signal is a frequency-hopped, minimum-shift keyed (FH/MSK) sig- 
nal. An FFT of the signal 1s shown in Figure | to illustrate the frequency-domain repre- 
sentation of the signal. The shape of the power spectral density function is approximated 
by the squared magnitude of the FFT. To simplify the algorithm, the magnitude of the 
FFT is used to represent frequency-domain signals in this thesis. The signal hops over 
5.8 MHz of spectrum, from 6.9 MHz to 12.7 MHz. Table 1 shows the approximate fre- 


quency and order of each of the seven frequency hops in the signal. 


Approximate | Hop 
Carrier Number 
Frequency 
of Hop 


[2.07 Miz [6 


Table 1. Frequencies and sequence earl: of re in FH/MSK signal. 
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Figure 1. Fourier transform of FH/MSK transmitted signal without interference. 


B. INTERFERENCE WAVEFORMS 
The interference signal is a composite waveform consisting of two separate sig- 
nals. The first file contains a continuous wave (CW) signal that has defined spectral 
components at 21 frequencies in the operating bandwidth of the frequency-hopped signal. 
This signal produces many spectral spikes, which closely resemble a frequency-hopped 
signal in the frequency domain. A plot of the magnitude of an FFT of the CW signal is 
4 


shown in Figure 2. The second signal 1s a binary-shift keyed signal (BPSK) with a car- 
rier frequency of 10 MHz. Figure 3 contains a plot of the power spectral density of the 
BPSK signal. Summing the signals together in the time domain created the total interfer- 


ence signal. A plot of the power spectral density of the total interference signal 1s shown 


in Figure 4. 
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Figure 2. Fourier transform of continuous wave signal. 
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Figure 3. Fourier transform of BPSK signal. 
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Figure 4. Fourier transform of total interference signal. 
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When all three signals are transmitted over the same time period, the spectrum 
becomes crowded and the frequency hops cannot be resolved directly. As shown in Fig- 
ure 5, the frequency hops are buried amongst the BPSK interference or are otherwise dif- 
ficult to distinguish from the CW signal. From Figure 5 it is evident that a filter-based 
detection method, like those described 1n [1] and [2] and shown in Figure 6, would be 
overwhelmed by the interference in the band. Each filter is meant to detect energy at a 
specific frequency, but interference signals that span the entire bandwidth of the desired 
signal and have a large amount of energy will cause each branch of the detector to indi- 


cate a frequency hop present at that frequency. 
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Figure 5. Frequency domain plot of total received signal. 


This 1s a significant problem because the frequency hops are essentially being 
jammed by other modulated waveforms without flat power spectral densities. From [5], 
additive white Gaussian noise is assumed to have a flat power density across all frequen- 
cies in the band. Gaussian noise that is uniform across all frequencies in the band 1s easy 


to estimate mathematically for any frequency hop in the spectrum. The spectral varia- 


tions of a communications signal such as the interference waveform mentioned in this 


thesis cannot be approximated as Gaussian noise. 






 Guadratic 
Detector 





a Quadratic: 
| ae 


Figure 6. Channelized receiver followed by likelihood-ratio test (From Ref. 1.). 


The signals used in the simulations were described in this chapter. From an ex- 
amination of past works, to the best of the author’s knowledge this is the first time that 
frequency-hopped signal detection has been studied in an environment containing non- 
Gaussian noise. A FH/MSK signal is the desired transmitted signal, but at the receiver 
the frequency-hopped signal is eclipsed by interfering BPSK and CW waveforms. A 
composite of the BPSK, CW, and FH/MSK signals is treated as the total received signal. 
In the next chapter, the use of the exponential averaging algorithm to form an estimate of 


the interference signals is described. 


Hl. INTERFERENCE ESTIMATION 


In this chapter efforts to estimate the interference present in the operating band 
are discussed. An exponential average of FFT windows was used to create an approx1- 
mate frequency-domain representation of the interference signal. The parameters of the 
exponential average, such as weight factor, FFT window size, and the number of win- 
dows in the average, are discussed. 

A. EXPONENTIAL AVERAGE 

The term interference estimation refers to a frequency-domain representation of 
the interference signal. In order to recover the frequency-hopped signal embedded in the 
interference signal, an acceptable estimate of the composite BPSK and CW signal is re- 
quired. As a first step, the components of the frequency-hopped signal must be removed 
from the total received signal spectrum. Using an exponential average of FFT windows, 
we can obtain an approximate picture of the interference spectrum. This algorithm 1s 


useful for averaging a discrete signal that has an unknown length. 


The interference estimate is created by exponentially averaging many small FFTs 
over time. The received signal, composed of the FH/MSK signal, as well as the BPSK 
and CW interference waveforms, was broken into N-point windows. The FFT of each 
window was computed and the magnitude of the FFT was taken. (When an FFT window 
is referenced 1n this thesis it can be assumed that it is the magnitude of a fast Fourier 
transform being described.) The FFT windows were used to form an estimate according 
to an exponential average. If the number of frames was known before averaging the win- 
dows, it would be simple to apply an equal-weight average to each window of 1/L, where 
Lis the number of windows used to form an estimate of the interference. The length of 
the signal is not known, however, and thus an exponential average is used to form the 1n- 


terference estimate. The averaging algorithm is represented mathematically by 
X[i]=(1-R)- X[i-1]+ R- Cli] (1) 

where X [7] 1s the i-th averaged window in the signal known as the reference frame, R 1s a 

weight factor between 0 and 1, Cli] 1s an FFT of the current window in the signal, 


i=1,2,K ,L, andL =1,032,192/N. The reference is multiplied by a factor of 1— R and 


9 


added to the current frame, which is weighted by R. The sum then becomes the new ref- 
erence frame. From (1), it can be shown that three important variables affect the final 
estimate: the FFT window size, N; the number of FFTs used to create the estimate, L; and 
the weight factor, R. 
B. FFT SIZE 

The size of the FFT affects not only the number of frames used in the estimate but 
also the resolution of the estimate. Since the final estimate will be the same size as the 
windows used to form the estimate, choosing the correct N 1s a key step in forming an 
interference estimate. As N increases, the resolution improves but the number of frames 
used to form the estimate, L, decreases because the signal is of finite length. A small N 
produces a well-averaged estimate but does not have enough detail to recover the fre- 
quency-hopped signal. For instance, a window size of 4,096 points lacks the resolution 
to accurately represent the detail in the spectrum. The resolution for the plot is 


Np el 0071 (2) 
NT. N_ 4,096 


where N = 4,096. An FFT of the received signal with only 4,096 points is shown in Fig- 
ure 7 to illustrate the difficulty in using an estimate with considerable error. The FFT has 
the general shape of the interference signal but contains many small nulls as a result of 


having only 4,096 points to represent the spectrum. 


Since 4,096 1s a factor of the signal length, 1,032,192, the FFT window size was 
increased in increments of 4,096 until satisfactory resolution was achieved. Integer mul- 
tiples of 4,096 ensure that the entire signal 1s averaged. If the window size does not di- 
vide evenly into the length of the signal, then a portion of the final samples in the signal 
may be excluded from the average. Thus, it is important to use an integer number of 


frames in the average. 
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Figure 7. A 4,096-point FFT of the received signal. 


The signal file is a finite length and, as a result, as the FFT window size grows, 
the number of FFTs that are used to calculate the average decreases. This results in an 
interference estimate that is formed from a small number of FFTs. Examination of win- 
dow sizes from 8,192 points to 147,456 points proved that more than one FFT size is 
suitable to produce a detailed estimate of the interference. Any size FFT window can be 
used to form an estimate of the interference as long as it meets the requirements set by 

ST (3) 


res 


where F; is the sample rate, Tyo, 1s the duration of a frequency hop, Af... is the resolution 


in the plot, and N is the FFT window size. This limit ensures that the window size N does 
not exceed the number of samples taken in a single frequency hop period. Averaging a 
smaller number of frames than there are frequency hops 1n the signal gives an invalid es- 
timate by not allowing enough frames to average out the frequency hops. Moreover, 


computing the FFT of a large frame size can be computationally intensive and limit the 


1] 


speed at which an estimate can be formed. In Table 2, a list of window sizes is compared 
with the number of frames used to form the estimate and the approximate resolution. 


Number Of |Approximate 
Windows Resolution, 


Averaged Afres [HZ] 
4,096 12,207 
8,192 6,104 





63 3,502 





Table 2. | Comparison of the relative FFT size, the number of FFTs averaged together 
to form the estimate, and the approximate resolution. 





For this thesis an FFT window size of 24,576 points was chosen. It should be 
noted that although this window size was not chosen arbitrarily, other window sizes 
would have sufficed. This frame length provided enough detail to depict frequency hops 
but also allowed the estimate to be calculated from 42 frames. As the number of frames 


used in the average increases, the frequency-hopped signal is diminished more. 


To ensure this window size meets the requirements set forth in (3), further inspec- 
tion is required. Each signal was sampled at 50 MHz. The sampling period then 1s 


1 1 
T =—=——_ =20ns. 4 
‘ F 50-10° “) 


From the sampling period and the frame length of 24,576 points, we can calculate the 
window duration with 
(24,576 samp/s)(20 ns) = 491.52 ys. (5) 


With seven hops in the time signal and the sampling rate of 50 MHz, the length of the 


time signal is 


(1,032,192 samp)(20 ns/samp) = 20.64384 ms. (6) 
The duration of a hop is then 
20.64384 ms/7 =2.94912 ms/hop, (7) 
making the hop rate 
1/2.94912-10° s » 339 hops/s. (8) 
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From (5) and (7), a window size of 24,576 points is well under the limit set in (3) since 
491.52 us < 2.94912 ms. (9) 

With a FFT size selected, the weight factor variable in the exponential average, R, can be 
selected. 
C. WEIGHT DETERMINATION 

The weight factor is the most important attribute of the exponential averaging al- 
gorithm. There is less flexibility in choosing a weight than in deciding on an FFT win- 
dow size. A change in the weight factor, as little as 0.05, can provide noticeably different 
results. From (1), it can be seen that the weight factor is the parameter that emphasizes a 
particular frame, thus regulating the contribution a single frame has on the overall esti- 
mate. A small weight factor accentuates the early-time hops in the sequence because 
each new frame only contributes a small amount to the estimate; hence, the effect of the 
early-time hops does not rapidly diminish. As a consequence of a small weight, a frame 
containing a frequency hop receives a large immediate reduction in power when multi- 
plied by R, but the residual power decays more slowly than when a larger weight is used. 
A large weight factor emphasizes the late-time hops, causing their immediate contribu- 
tion to the estimate to be large. Although the frames containing the late-time hops are 
accentuated with a large weight, the residual power of the frames containing early-time 
hops is quickly attenuated by the factor of 1— R. Fora large R, thel—R factor signifi- 
cantly deemphasizes the lingering frequency hop power left in the reference frame of the 
interference estimate. A balance must be struck between these two conditions by select- 
ing a weight factor that allows the contributions from individual hops to be as uniform as 


possible. 


Figures 8, 9, 10, and 11 are plots of the interference estimates for different 
weights. Each weight has tradeoffs between slightly degrading many of the frequency 
hops from the estimate and reducing hop power for all but a few hops. For instance, in 
Figure 8 the effects of using a weight factor of 0.05 to form the estimate can be seen. 
This weight creates an estimate with a nice shape of the interference spectrum but does 
not average out a single frequency hop. All seven hops can be seen in the estimate, and 
those at 12.07 MHz and 6.9 MHz still have significant amounts of power. Although the 
hop power 1s noticeably reduced, this weight 1s not successful in removing the frequency- 
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hopped signal from the spectrum. With significant portions of the desired signal remain- 


ing in the estimate, there is a risk that none of the hops will be recovered. 


A weight factor of 0.15 was used to form the estimate shown in Figure 9. This es- 
timate has more encouraging results than an estimate using a weight of 0.05. The contri- 
bution of the frequency hop at 9.58 MHz is more significant because it 1s the last hop in 
the sequence and thus weighted by R the fewest number of times. The tradeoff for having 
a large hop at 9.58 MHZ in the estimate is a greater reduction in the power of other fre- 
quency hops. The hop at 12.07 MHz 1s large in this estimate, but little can be done to re- 


duce that hop without severely impacting others because it is in the middle of a spectral 


null. 
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Figure 8. Interference estimate formed from weight of 0.05. 
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Figure 9. Interference estimate formed from weight of 0.15. 


A weight factor of 0.25 provides results similar to that of a weight of 0.15. The 
last hop in the sequence has a larger power in the estimate, but the frequency hop at 6.9 
MHz is reduced below the power level of the BPSK signal. A plot of the power spectral 


density of the estimate is shown in Figure 10. 


A plot of the estimate using a weight of 0.5 is shown 1n Figure 11. This weight is 
not a good number to use in the exponential average. It forms half of the estimate from 
the last frame and the other half of the estimate from the previous 83 frames. This is un- 
desirable because the final hop at 9.58 MHz has approximately 10 dB more power than 
the BPSK signal at that frequency, making it a considerable part of the interference esti- 
mate. When the interference estimate is used to recover the frequency-hopped signal, the 
last hop will be treated as interference and will not be recovered. Furthermore, the esti- 
mate provides a rougher estimate of the BPSK spectrum than estimates obtained with a 
smaller weight. The variation in the shape of the BPSK spectrum is greater, resulting in 


an estimate that 1s not as smooth as the estimates obtained with smaller weights. This 
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deviation could provide problems when this estimate is used to resolve the frequency- 
hopped signal. Although the estimate 1s not perfect, a reduction in the power of the fre- 
quency hops to as close as possible to the peak power spectral density of the BPSK signal 
is desired. A weight of 0.5 reduces five hops to that level but is unsatisfactory because 
the estimate places too much emphasis on the last hop and the shape of the BPSK spec- 


trum 1s noticeably coarse. 


In Table 3, a list of weights and the respective number of hops left in an estimate 
formed with that weight 1s shown. Smaller weights leave the most frequency hops 1n the 
estimate, and the number decreases as the weight grows larger. The table leaves out the 
qualitative effects of each weight factor though, not showing that, as weights increase, the 
powers of the remaining hops in the estimate increase. The hops left after using a weight 
of 0.5 have a large amount of power in the estimate, and the powers of residual hops only 


grow when larger weights, such as 0.7, are used. 


T 


Magnitude of FFT 





10° 
5 6 f a 9 10 11 12 13 14 15 
Frequency (MHz) 
Figure 10. Interference estimate formed from weight of 0.25. 
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Figure 11. Interference estimate formed from weight of 0.5. 


Number 
Of Hops 
Remaining 





Table 3. © Number of frequency hops remaining in interference estimate for given 


weight. 


From the plots of the interference estimates, a weight between 0.15 and 0.25 was 
determined to be the ideal factor to use in the exponential average algorithm. Weights in 


this range stress individual frames as equally as can be expected for an algorithm that 
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does not know the duration of the signal at the start of executing the algorithm. To de- 
termine the best weight factor for recovering the frequency-hopped signal, the different 


estimates must be tested to determine which estimate recovers the most frequency hops. 


In this chapter, the development of a frequency-domain estimate of an interfer- 
ence signal eclipsing a desired frequency-hopped signal was discussed. With an expo- 
nential average of FFT windows, a frequency-hopped signal can be removed quite well 
from the estimate. With a weight factor of 0.25 and a window size of 24,576 points, an 
accurate estimate of the interference signal spectrum was obtained. In the next chapter, 
we examine various methods that utilize the interference estimate for detecting the fre- 


quency-hopped signal embedded in the interference. 
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IV. FREQUENCY HOP RECOVERY 


With acceptable estimates of the interference signals, a method was developed for 
resolving the frequency-hopped signal buried within the noise. Two approaches were 
examined: spectral subtraction and spectral division. The goal of the two techniques 1s to 
exploit the difference between the interference estimate and the received signal. Ideally, 
that difference should be only the desired frequency-hopped signal. Since real signals are 
not easy to estimate, the objective was to find the best case from either of the two meth- 
ods. In order to accomplish the mathematical calculations of array values, the power 
spectral density of the received signal had to be the same length as the estimate. Once the 
estimate was formed, a frame-by-frame calculation was made. The results were then 
summed to create a final spectral picture that was 24,576 points long. 

A. SPECTRAL SUBTRACTION 

Spectral subtraction assumes that the frequency-hopped signal has more power 
than the interference signal at each of the hop frequencies. The subtraction algorithm op- 
erates on the principle that a hop with 1-2 dB more power would be recognizable once 
the interference was removed. Spectral subtraction, however, turned out to be an unsuc- 


cessful method for recovering frequency hops. 


Since the estimate is imperfect, it 1s not able to calculate the actual power levels at 
each frequency; only the general shape of the spectrum 1s preserved with the exponential 
average. As aresult, spectral subtraction highlights the difference between assumed and 
actual power levels of the interference signal. The result, known as the difference spec- 
trum, then retains the same general shape as the interference signal. No separation be- 
tween the frequency-hopped signal and the interference signal was obtained, and the 
FH/MSK signal remained embedded 1n the interference signal. In Figure 12 a plot of the 
difference spectrum is shown. It 1s noteworthy that at 9.58 MHz and 12.07 MHz hops are 
recognizable but all other hops could easily be confused with noise and cannot be accu- 


rately resolved. 
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Figure 12. Power spectral density of spectral subtraction signal, weight fac- 
tor R = 0.25. 


The failure of spectral subtraction to resolve the frequency-hopped signal high- 
lights the fact that numerical power levels in the estimate are not the important factor. 
The comparative shape between the estimate and the received signal are the keys to de- 
tecting the frequency-hopped signal. Normalization through division exploits the shape 


of the estimate rather than comparing the power levels in the two spectrums. 


B. SPECTRAL DIVISION 
As a byproduct of the spectral subtraction results, the idea of division was devel- 
oped. If the frequency-hopped signal can viewed relative to a noise floor and not the in- 


terference signal as with the original composite received signal, then the power difference 


at the hop frequencies should be noticeable. 


A frame-by-frame division of the received signal by the estimate was computed. 
Again the results were summed, and a 24,576-point power spectral density showed the 


difference in shape of the received signal and the estimate rather than the difference in 
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power. The frequency-hopped signal had then been normalized against the interference 
signal and was able to be viewed as if it were the only signal in the band when received. 
A type of receiver that implements division of the received signal with the noise power is 
known as a noise-normalized receiver. Considerable analysis has been done in the past 
on these types of receivers, and the calculated theoretical performance asserts that noise- 
normalized receivers significantly reduce bit error rates. Previous work on noise- 
normalized receivers in [6], [7], [8], and [9] assumes that the received noise is uniform 
and Gaussian, a considerably different case compared to the interference signals consid- 
ered in this thesis. None of these references discusses implementation of the receiver in 
hardware or simulation; all analysis is theoretical. The estimate used for normalization in 
this thesis is formed over the length of the signal, while in [7] and [9] the noise estimate 
is formed from interference present only during the hop period. With the estimate de- 
scribed in Chapter 3, a receiver is able to normalize the desired signal even with non- 


uniform noise power. 


The results of the spectral division algorithm were considerably better than those 
achieved with subtraction. Figure 13 shows the normalized spectrum. From this spectral 
picture it 1s quite evident that frequency hops are present in the received signal. Although 
not all the hops are present, division with the estimate yields a majority of the frequency 
hops. Using a weight of 0.25, five of the seven transmitted frequency hops are clearly 
defined in the spectrum. The two missing hops are the sixth and seventh 1n the sequence. 
The hops are numbered with the order in which they were transmitted. Hop 6 1s buried in 


the noise floor and hop 7 appears to be negative relative to the noise floor. 


The estimate formed with a weight factor of 0.15 was also used to normalize the 
received signal. The smaller weight did not yield as many frequency hops as the estimate 
that used 0.25. Hops 6 and 7 are unable to be recovered, but the fifth hop in the sequence 
is too close to the noise floor to be resolved as well. In Figure 14 the normalized spec- 
trum is shown after 0.15 was used for the estimate weight factor. The noise floor in that 
spectrum has less power than the floor in Figure 13 for a weight of 0.25. Although this is 


an attractive quality, 0.25 still allows the recovery of the most frequency hops. 
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Figure 13. A plot of the power spectral density of the normalized spectrum, R = 0.25. 
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Figure 14. A plot of the power spectral density of the normalized spectrum, R = 0.15. 
pis 


In order to confirm that estimates formed using other weight factors would not 
perform as well as a weight of 0.25, other weights such as 0.5 and 0.05 were tested. As 
expected, these weights did not provide an estimate accurate enough to recover the fre- 
quency-hopped signal as well as a weight of 0.25. In Figure 15 a plot with the normal- 


ized spectrum for R =0.05 is shown. Although the spectrum is very clean, only four 


hops were recovered. 
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Figure 15. A plot of the power spectral density of the normalized spectrum, R = 0.05. 


Increasing the weight does not cause more frequency hops to be recovered either. 
As shown in Figure 16, with a weight of 0.5, six hops can be seen but are difficult to dis- 
tinguish from the noise floor in the spectrum. Both hops five and six could be mistaken 
for noise spikes, and it would be difficult to recover the demodulated signal amidst the 
higher noise levels. As the weight factor increases, the noise variance relative to the fre- 
quency-hopped signal increases and eventually overcomes the reduction in noise power 


provided by normalization. A weight factor of 0.5, as well as larger weights, caused the 
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noise variance to increase because much of the estimate is formed from the last window 
in the received signal. With so much emphasis placed on the last frame, the estimate 1s 
less of an average and instead relies on the final window being similar to every other 
window in the received signal. It is not, however, and much of the noise arises because 


the average has not been able to smooth out the shape of the estimate. 
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Figure 16. A plot of the power spectral density of the normalized spectrum, R = 0.5. 


From the previous figures, it is evident that a weight of 0.25 resolves the most 
frequency hops without exciting the noise floor too much. This makes recovery and de- 
modulation of a frequency-hopped signal an easier task. Using this weight, we assumed 
that all but the last two frequency hops transmitted could be recovered using the noise- 
normalization technique described. In order to determine how this technique performs 
with longer received signals, a transmitted signal was used that had a noise file appended 


after the transmission of the frequency-hopped signal had ended. This is discussed in the 


next section. 
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C. NOISE-NORMALIZATION WITH EXTENDED RECEIVED SIGNAL 

Under the theory that forming the estimate from a longer signal contributes to the 
detection of late-time frequency hops, a simulation was run to determine if all seven fre- 
quency hops in the signal could be recovered. A noise file was appended to the end of 
the received signal to give the estimate more time to normalize the received signal. 

1. Extended Signal with Appended Noise 

A received signal file was formed from the original received signal from Chapter 
II, and then a 1,032,192-point file containing only the interference waveforms was ap- 
pended at the end. This file was then twice the size of the original received signal. This 
signal simulates a frequency-hopped waveform being received amidst interference wave- 
forms but, when transmission of the frequency-hopped signal stops, the receiver still con- 
tinues to receive the interference waveforms. A spectrum of this signal looks the same as 
the spectrum shown in Figure 4. The time signal 1s different, however, than the original 


received signal. A plot of the received signal in the time domain is shown in Figure 17. 


FH/MSK signal + interference signal 


Interference signal only 
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Figure 17. Time domain plot of extended received signal. 
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2. Simulation results of extended received signal 

After normalization, the spectrum of the extended received signal (Figure 18) 
looked very clean. All seven frequency hops are clearly shown. From this simulation, it 
is evident that, as a discrete signal grows in length and L increases, frequency hops re- 


ceived in the i-th window of the signal, where 7 is small compared to L, can be resolved 


with more clarity. 
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Figure 18. Plot of power spectral density of extended received signal after normaliza- 


tion. 


As aresult of the increase in L, the contribution from earlier windows is reduced 
in the estimate. The window that contains the seventh hop in the sequence is no longer 
the last frame, and thus no longer a large portion of the estimate. This was the reason the 
seventh hop could not be resolved in previous simulations. A window that encompasses 
the last frequency hop, the 42"° of 42 windows, is now the 42™ of 84 windows. The 42" 
window had previously been multiplied by R=0.25. The weight factor on that individ- 


ual frame is now 
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(R)(1-R)”) =(0.25)(0.75)*? ~1.4-10°, (10) 
The contribution of the 42" frame to the interference estimate is virtually zero, thus cre- 
ating an estimate with no frequency hop components. From these results, it is clear that 
the algorithm should be able to resolve the first k -2 hops inak hop signal. This 1s ex- 


amined in the next section. 


D. NOISE-NORMALIZATION WITH EXTENDED TRANSMITTED 
SIGNAL CONTAINING FOURTEEN FREQUENCY HOPS 


Although noise-normalization was able to resolve most of the frequency hops in 
the received signal, another simulation was devised to test ifk —2 hops in a k-hop signal 
could be detected. A signal with fourteen frequency hops was created for this simulation. 

1. Signal with Fourteen Frequency Hops 

Since no signal with a greater number of frequency hops was available, a file was 
created using the original frequency-hopped signal. The seven hops were shifted in fre- 
quency to new, unique carriers and concatenated in time with the original FH/MSK signal 
to provide a signal with k =14 frequency hops. To achieve this, a frequency shift was 


applied to create seven new hop frequencies. Using the Fourier transform property [10] 


ell! <= 5(f — fy) (11) 


J2m fot 


to shift the frequencies, the time domain signal was multiplied by e , where 


fy =1.3 MHz. A frequency-domain plot comparing the original seven hops and the 


seven shifted hops is shown in Figure 19. 


The spectrum of the complete fourteen-hop signal is shown in Figure 20. The hop 
sequence for the signal is shown 1n Table 4. Because the fourteen-hop signal was com- 
prised of two concatenated signals, the first seven frequency hops are the same as the 
original signal and the last seven hops are transmitted in the same sequence as the first 


seven. The concatenated signal is 2,064,384 points long. 
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Figure 19. Plots comparing spectrums of original and shifted frequency-hopped sig- 
nals. 
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Figure 20. Power spectral density of frequency-hopped signal with fourteen hops. 
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Carrier Hop 
Frequency | Number 
[MHz] 





Table 4. | Frequencies and sequence order of hops in extended, 14-hop FH/MSK signal. 


2 Extended Signal Detection Results 

As the weight increases, more noise is left in the normalized spectrum. Increased 
weight factors also place greater emphasis on later hops in the sequence, accentuating 
hops that would not be detected with smaller weights. A weight of 0.3, slightly larger 
than the weight factor of 0.25 used in Section B, provided the best normalized spectrum. 
Although a weight factor of 0.25 was satisfactory, the hop at 5.6 MHz was clearer when 
R=0.3. For every seven hops transmitted, six were resolved. In this simulation twelve 
of the fourteen hops that were transmitted are visible in the normalized spectrum. A plot 


of the power spectral density of the normalized signal is shown in Figure 21. 


Similar to the results discussed in Section B, small weights such as 0.05 and 0.15 
do not resolve as many frequency hops as larger weight factors. Any weight factor above 
0.3 in this simulation made it too difficult for certain hops to be differentiated from the 


noise floor. 


From these results it 1s clear that normalization of the spectrum allows six of 
every seven hops to be resolved if a weight of approximately 0.25 to 0.3 is used. This 
upper limit of six hops per every seven 1s a positive result when compared to interference 
waveforms remaining in the received spectrum and no hops being detected. It is impor- 


tant to note that all of the first seven hops in the sequence were resolved. The only two 
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hops not resolved were numbers thirteen and fourteen in the sequence. From this it can 
be assumed that the earlier a hop is in the sequence, the more likely that it will be re- 
solved. This follows intuitively because only the late-time hops were not averaged out in 
the estimate. Since hops thirteen and fourteen are significant in the estimate, they do not 


stand out when the received signal is normalized by the estimate. 
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Figure 21. Normalized spectrum of extended received signal, R = 0.3. 


In this chapter, several simulations were used to test methods for resolving a fre- 
quency-hopped signal. Subtraction of the estimate from the received signal proved to be 
an inadequate method for detecting the frequency-hopped signal because too much noise 
remained in the difference spectrum. Spectral division, or noise-normalization, however, 
is a technique that successfully resolved the frequency-hopped signal embedded 1n inter- 
ference. With this technique, k —2 hops in a k-hop signal can be recovered. The next 


chapter focuses on the performance of noise-normalization with regards to the bit error 


rate and interference power reduction. 
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VV. BIT ERROR RATE SIMULATION 


A simulation measuring the effectiveness of the noise reduction algorithm was 
tested with a new transmitted signal. The goal was to determine the reduction in interfer- 
ence signal power necessary to achieve the same probabilities of bit error (1.e., bit error 
rates (BERs)) for a received signal, viewed after the channel noise had been normalized. 
A. SIMULATION DESCRIPTION 

Three simulations were needed in order to characterize the effects of the normali- 
zation technique. A frequency-hopped, binary frequency-shift keyed signal (FH/BFSK) 
was transmitted in a channel with various types of noise. If perfect de-hopping 1s as- 
sumed at the receiver, then the received signal 1s a BFSK signal embedded 1n interfer- 
ence. To test this, additive white Gaussian noise (AWGN) was added to the channel as 
well. The first simulation measured the BER fora MATLAB-generated BFSK signal in 
AWGN. Figure 22 shows a block diagram of the simulation. The second simulation was 
of the BFSK signal transmitted through a channel that contained BPSK and CW wave- 
forms as well as AWGN. The same BPSK and CW waveforms described in Chapter II 
were used for the simulation. The block diagram shown in Figure 23 describes this simu- 
lation. The third simulation measured the BER for a BFSK signal transmitted over a 
channel with AWGN and a normalized interference signal. This interference signal was 
created by computing the inverse fast Fourier transform (IFFT) of a normalized compos- 
ite BSPK and CW interference signal. This signal is known as the normalized interfer- 
ence signal. A block diagram describing the simulation is shown in Figure 24. In this 
simulation the performance of a communications signal sent through a channel with nor- 
malized interference power levels is examined. With a normalized interference signal, 
noise power levels should be significantly lower than the simulation shown in Figure 23. 
The reduction in interference power was measured by comparing the bit error rates from 


the simulations before normalization and after normalization. 
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Figure 22. Block diagram of simulation measuring bit error rate for received signal in 
channel with AWGN. 
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Figure 23. Block diagram of simulation measuring bit error rate for received signal in 
channel with AWGN and interference signal. 
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Figure 24. Block diagram of simulation measuring bit error rate for received signal in 


channel with AWGN and normalized interference signal. 
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B. SIMULATION RESULTS 

The signal transmitted through a channel with just AWGN was received with an 
approximate probability of bit error of 0.001. To determine the BER for the 65,536 
transmitted bits, a BFSK-receiver simulation, standard as part of the MATLAB commu- 
nications toolbox, mapped the received symbols to bits and compared each with the trans- 
mitted bits. With interference waveforms present 1n the channel, the BER degraded to 
approximately 0.44. This level of bit error rate is essentially receiver failure because the 
receiver is forced to guess, correctly choosing the transmitted symbol only half of the 


time. 


The probability of bit error for the channel with a normalized interference signal 
and AWGN was 0.001. Since this is the same as with only AWGN, the capabilities of 
the normalization scheme are best characterized by the amount that the algorithm reduces 
the power of the interference signal. In order to quantify this reduction, the power of the 
transmitted interference waveforms was reduced until the same bit error rate as the simu- 
lation with a normalized interference signal was achieved. The plot in Figure 25 shows 
how the bit error rate changed with a reduction 1n interference signal power. As seen in 
Figure 25, the bit error rate drops dramatically as the interference signal power 1s re- 
duced. In order to achieve the same bit error rates in the simulation as with normalized 
interference and AWGN, the interference signal power had to be reduced by 50 dB. It 
can be assumed that the noise-normalization technique could achieve the same reduction 


in interference signal power when used for a real-time reception. 


This bit error rate simulation provides a quantitative characterization of the noise- 
normalization effects. Use of the normalization technique yields bit error rates equal to 
those of a 50-dB reduction 1n interference signal power. From this, it can be assumed 
that noise-normalization reduces interference by 50 dB 1n the frequency-hopped band. 
Chapter VI uses this result to draw conclusions about the detection of frequency-hopped 


signals embedded in interference waveforms. 
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* BER of AWGN + normalized interference simulation 
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Figure 25. Plot of reduction in interference signal power vs. bit error rate along with 


plot of BER for simulation with noise signal containing AWGN and nor- 
malized interference signals. 
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VI. CONCLUSIONS 


A. CONCLUSIONS 

In this thesis a technique for resolving a frequency-hopped signal embedded in 1n- 
terference waveforms was described. The algorithm involves creating an estimate of the 
interference spectrum by using an exponential average of FFT windows to create an ap- 
proximate spectral picture of the interference waveforms. The variables involved in av- 
eraging the windows include the FFT size, the weight applied to each frame, and the 
number of frames used in the estimate. For a 24,576-point window, a weight between 
0.15 and 0.25 was determined to produce the best estimate of the interference spectrum. 
Two methods for resolving the frequency-hopped signal were explored as well. Division 
of the received signal by the estimate 1n the frequency domain presents a cleaner power 
spectral density of the frequency-hopped signal than does spectral subtraction. Fora 
weight of 0.25, five of the seven frequency hops in the desired signal were detected after 
normalization of the received signal. Through other simulations, it was determined that 
k —2 frequency hops of a k-hop signal can be resolved using this noise-normalization 
technique. Through bit error rate simulations, the normalization technique was deter- 
mined to have reduced the interference waveform power by 50 dB, having achieved the 


same bit error rates as a 50-dB reduction in noise power. 


This technique may be useful for systems implementing frequency-hopped com- 
munications such as military radios and networks as well as commercial frequency- 
hopped systems such as Bluetooth [11]. Although these systems do not always operate in 
environments that contain other waveforms which would act as jamming signals, the pos- 
sibility exists that they might. 

B. FUTURE WORK 

The technique presented has future research potential and should be explored on a 
deeper level. While the method was developed for CW and BPSK-modulated interfer- 
ence signals, other types of waveforms that might interfere with frequency-hopped com- 
munications should be explored. These interference signals include M-ary phase-shift 


keying (M-PSK), M-ary frequency-shift keying (W-FSK), quadrature amplitude modula- 
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tion (QAM), and orthogonal frequency-division multiplexing (OFDM). Testing in differ- 


ent operational environments may allow for more uses of normalization. 


Furthermore, implementation of the technique in a real-time environment using 
simulations or hardware is another area to explore. Such an implementation would make 
the technique more applicable to real-time communications systems. Actual demodula- 
tion and recovery of encoded messages, as well as performance evaluations in fading 


channels, are steps that could be taken in future work. 


More can be explored regarding the capabilities of noise-normalization by inves- 
tigating the correlation between FFT window size and weight factor and the quality of 


normalized signal each combination produces. 
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